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A fundamental constant in systems of unitary fermions is the so-called Bertsch parameter, the 
ratio of the ground state energy for spin paired unitary fermions to that for free fermions at the 
same density. I discuss how we computed this parameter as well as the pairing gap using a 
recently developed lattice construction for unitary fermions, by measuring correlation functions 
for up to 38 fermions in a finite box. Our calculation illustrates interesting issues facing the study 
of many-body states on the lattice, which may eventually be confronted in QCD calculations as 
well. 
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1. Introduction 

Fermions at unitarity is an idealized system which has an infinite 5-wave scattering length and 
zero interaction range or equivalently the two-particle s-wave scattering phase shift 8q = n/2. In 
the unitarity limit, the only relevant scale is the density (p) and the ground state energy is related 
to that of the non-interacting fermions by 



where the dimensionless parameter E, is called the Bertsch parameter. E, is of our particular interest 
because it is the unique parameter which relates zero temperature thermodynamic quantities be- 
tween unitary and free Fermi gas. Another important physical quantity is the pairing gap (A), the 
difference between the chemical potential and the minimum energy required to add one fermion 
to the unpolarized unitary Fermi gas. This quantity plays an important role in characterizing the 
superfluidity of the system. 

In the unitary regime the interaction is strongly attractive and non-perturbative calculations 
are required. Since the details of the microscopic interaction are unimportant, such a system can 
be effectively described by four-Fermi contact operators with two-component fermions in the La- 
grangian. Based on this theory, we recently developed a highly improved lattice calculation method 
to explore a large number of strongly interacting non-relativistic fermions [p. We briefly review 
this method in Sec. 2. There has also been a substantial number of numerical studies of the unitary 
fermions from the microscopic theory using the Quantum Monte Carlo [Q, H, S] and various other 
Monte Carlo techniques [||, 

The underlying idea of the lattice theory is that the correlation function decays with the ground 
state energy in the long imaginary time limit. In principle, this is true as long as the interpolating 
operator has some overlap with the ground state and the error propagation is moderate with time. 
However, practically it is very difficult to find the ground state in simulating many interacting 
particles. Because the overlap between the source and the ground state decreases exponentially with 
respect to the number of particles (AO and the accessible time range is relatively small for large N, 
such as our case, as well. One way of overcoming these issues is by constructing improved sources 
(or sinks), which have large overlap with the ground state and result in earlier plateaus in effective 
mass plots. In Sec. 3, we discuss our choice of sources in detail and demonstrate the improvement 
that can be achieved with an example. In Sec. 4, with these improvements we calculate the ground 
state energy of up to 38 unitary fermions in a finite box and report our preliminary results for E, and 
A. In our companion proceedings, we present a similar calculation for up to 20 unitary fermions in 
a harmonic trap [|8j]. 

2. Lattice construction for fermions at unitarity 

In this section, we briefly review the lattice construction for unitary fermions which is dis- 
cussed in Ref. ^ in detail. Starting with the lattice action for interacting non-relativistic two- 
component fermions which is first introduced in Ref. [Q], we construct the single fermion propa- 



Eunitary(p) — ^Efree(p), 



(1.1) 



gator 



K 1 (t; 0) = D 1 ( 1 - (t - 1 ) Vc)D~ 1 • • • D 1 ( 1 - (0) Vc)D 



(2.1) 
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where D 1 = 1 - V 2 /2M, 0(t) is a Z2 auxiliary field on each time-link produced stochastically, 
and C is a generalized version of four-Fermi coupling involving an effective derivative operators 
with appropriate coefficients. With open boundary conditions in the time direction and time-like 
links of s, one can show that the fermion determinant is independent of the auxiliary field and 
the full numerical simulation is quenched. By averaging out the 0s, one can achieve a four-Fermi 



interaction. The A-fermion propagator can be constructed by taking direct products of Eq. ( |2.1[ ), 
anti-symmetrizing, and then ensemble averaging. 

In order to minimize the lattice artifacts, we consider a couple of lattice improvements. First, 
we use D in momentum space by 

( e P 2 /(2M) S <A 

D(p)= ™' P l A . (2.2) 

I 00 , p > A 

where A = 71 is a hard momentum cutoff and p = |p|. With this definition of D, one can elimi- 
nate the discretization errors appearing in free fermion propagator. Second, we minimize the dis- 
cretization errors associated with the contact interaction on the lattice when we tune the coupling 
constants to describe the unitary limit. The key idea of our tuning procedure is the following : we 
tune the coupling constants by matching the low-lying energy eigenvalues calculated from the two- 
particle transfer matrix & in our lattice theory to those of a finite volume and continuum theory 
for /?cot5o = using Luscher's formula. In principle, with these tuning methods we can eliminate 
both discretization errors and finite volume effects. Actually, in the unitary limit these two errors 
are equivalent since the only relevant dimensionless quantity is b/L, where b is the lattice spacing 
and L is the size of the lattice. 



3. Af-particle correlation function and two-particle source 

Since we are working in the canonical ensemble, we consider a two-component system of 
A| spin up and N\, spin down fermions. To satisfy Fermi-Dirac statistics, fermions of the same 
species have different quantum numbers. For the purpose of measuring correlators, we take the 
lowest momentum eigenstates of an uncorrected single free fermion as sources, such as pi = 
(0,0,0)2tt/L, P2 = (1,0, 0)271 /L, • • • in the Cartesian basis of a cubic lattice with size L. 

For the unpolarized fermi gas, the A-body correlation function CV(t) (A^ = N± = N/2) can be 
obtained from a Slater determinant of the N/2xN/2 matrix as 

Cjv(t) = (det(S(T))>, Sij(z) = (Y S ink\K- l (T;O)0K- l (r;O)\p iPj ), (3.1) 

where Xjfsink is an interpolating field of a spin up and a spin down fermion for the sink, and i,j = 
1,2, • --N/2. In the case of A^ = N± + 1, we replace the Mi row by (p^|*r _1 (T;0)|py) with j = 
1,2,--.,A72. 

The sinks may be chosen as the lowest momentum eigenstates of an uncorrected single free 
fermion as well, which corresponds to the ground state of the non-interacting fermions How- 
ever, the two-particle 5-wave solutions to the Schrodinger equations for unitary and non-interacting 
fermions behave like 1 jr and 1 , respectively, where r is the relative distance between the two 
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Figure 1: Effective mass plots for L — 8, M = 5, and 
N = 8 (iVj- =Ni = 4). The purple triangles are ob- 
tained using the lowest momentum eigenstates for 
a uncorrected single free fermion while red circles 
and green squares are obtained using correlated two- 
particle states for spin up and down as sinks. The 
black dashed line correspond to the ground state 
energy for non-interacting fermions. b = 0.5 and 
A = 1 .01 are used for red and green, respectively. 



Figure 2: Effective mass plots for L = 12, M = 5, 
and N =\0 (N^ = N± = 5). The red circles and blue 
triangles are results using the sinks in Eq. (3.2) with 
b = 0.5 and b = 1 .0, respectively. The black squares 
are obtained using the Multi-Correlation function 
Prony Method with the data from these two sinks, 
where the parameter Tj and Tw are introduced to im- 
prove stability Jl0|]. 




fermions. In order to have better overlap of the interpolating field with the ground state of N uni- 
tary fermions in a finite box, we choose the sink for one spin-up and one spin-down fermion as a 
two-particle state given by 



VW(P): i » " P + \ (3-2) 
p = 



where p is the momentum of each fermion in the center of momentum frame, p = |p|, and b and 
\j/Q are tunable parameters. Ignoring the p ^ component, Ysmkiv) Wlt h b = corresponds to the 
wave function of the unitary fermions in an infinite volume and continuum theory. The result is not 
sensitive to Yo and we consider Yo = 100 in this work. 

As an example of the wavefunction overlap problem in many particle calculations and the 
importance of using a proper choice for the interpolating operator, we plot the effective masses with 
three different sinks for L = 8, M = 5, and N = 8 in Fig. |l[ Around % = 1 8 — 25, the effective masses 
with two different correlated two-particle sinks, red (corresponding to the operator in Eq. (3.2)) and 
green (corresponding to the eigenstate of the two particle transfer matrix with a tunable parameter 
A), have clean plateaus with consistent fit values while the effective mass with the uncorrected 
single particle sink (purple) is much higher and has no clear plateau. This implies that the purple 
doesn't reach to the ground state, yet. For T > 26 all of the effective mass plots become very noisy. 
This result shows that it is very difficult or even impossible to see the unitary fermion ground state 
using sinks constructed from single particle momentum eigenstates. 
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Figure 3: (Preliminary) Bertsch parameter as a 
function of N for V = 8 3 (blue), 12 3 (black), and 
16 3 (purple). 



Figure 4: (Preliminary) Bertsch parameter as a 
function of N for V = 12 3 . The solid line represents 
the fitting range and fit result for % . 



4. Preliminary Results 

In this work, we calculate the ground state energy of N < 38 unitary fermions and extract pre- 
liminary results for the Bertsch parameter and pairing gap. As discussed in the previous section, the 
plateaus in the effective mass plots are much clearer when the sinks are constructed from correlated 
two-particle states. Since we have different sinks for the calculation, we can obtain slightly longer 
and better plateaus using the Multi-Correlation Prony Method [[k]] with two different sinks (Fig. ^|). 
After blocking the data, we perform a Jackknife resampling procedure and then a % correlated fit 
of the plateau obtained using the Prony Method. The statistical error is obtained by varying the 
fit parameter to give x 2 = Zmin + 1 while the fitting systematic error is obtained by varying the 
end points of the fitting range by ±2, ±1,0. The quoted errors represent only the combination of 
the fitting statistical and systematic errors in quadrature and does not include the systematic errors 
from finite volume effects or discretization errors. 

As discussed in Sec. 2, in principle we expect that there should be no volume dependence in 
the Bertsch parameter. However, since we only tune a finite number of the coupling constants, the 
deviation of p cot 8 from zero gets larger with p, resulting in the presence of discretization errors 
for small V and large N. In Fig. |5[ we calculate ^ for N < 14 unitary fermions in a periodic box 
with V = 8 3 , 12 3 , and 16 3 . For N > 14, with V = 8 3 is slightly greater than that of the other two 
volumes. This result shows that V = 8 is too small to calculate the correct t, for N > 14. 

The results for the Bertsch parameter for N < 38 and V = 12 3 are shown in Fig. ^f. We find 
that £, is not independent of N for the first shell while we don't see the N dependence in the second 
shell with given statistics. A fit to the second shell gives ^ = 0.412(4) (preliminary). The value of 
(?; has been obtained from various theoretical calculations and ultra-cold atomic experiments. For 



example, the recent QMC calculation gives E, = 0.40(1) [11] while one of the recent experimental 
measurements shows ^ = 0.39(2) and = 0.41(2) using two different extinction methods [[T2|]. 
We also perform the simulation for the slightly polarized unitary fermi gas, = N± + 1. As 



1 We also did the simulation for L = 16, but for the given computer resources we did not have sufficient statistics to 
calculate the Bertsch parameter for the second shell. 
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Figure 5: (Preliminary) Ground state energies for 
N < 38. The red circle and blue square represent the 
ground state energies for unpolarized and slightly 
polarized unitary fermions. 



Figure 6: (Preliminary) Pairing gap in unit of \lf ree 
for N < 37, where pLf ree is the free Fermi energy. 
The solid line represents the fitting range for the ex- 
tracted value of A. 



discussed in Ref. [0], the ground state energy for odd N depends on the momentum of the unpaired 
spin up fermion. For N < 37, the minimum energies are obtained by placing the unpaired fermion 
on the first momentum shell, and used for Fig. || and the gap calculation. In Fig. ||, we plot the 
ground state energy of both the unpolarized and the slightly polarized unitary fermi gas. The 
staggering between even and odd indicates a non-zero pairing gap. Using this data, we calculate 
the pairing gap as 

A = E (N-f,Ni) 1 2 V t (4.1) 

and the results are shown in Fig. ^. By fitting the data in the second momentum shell, we obtain 
A = 0.52(1). The most recent QMC calculation [Q] and ultra-cold atomic experiment [|13|] reported 
A = 0.50(3) and A = 0.43(3), respectively. 



5. Conclusion 

In these proceedings, we explored up to 38 unitary fermions in a finite box, and extracted 
preliminary values for the Bertsch parameter and pairing gap. The absence of volume dependence, 
a consequence of the tuning method for the four-Fermi contact interaction, allows us to study the 
infinite volume, continuum properties of unitary fermions with a moderately sized lattice. Our 
analysis also shows that a careful choice of the interpolating operators, which are designed to have 
large overlap with the state of interest, plays a crucial role in simulating the interacting many- 
particle system on the lattice. In the future, we hope to extend our simulations to larger numbers 
of unitary fermions using a recently developed tuning procedure 2 and more careful study of the 
correlated two-particle sinks. 

2 We have recently developed a new tuning method which has Hermitian, Galilean invariant and analytic interaction. 
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